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Q ■ Abstract 

O 

CN ■ I review the lattice simulations of light scalar mesons with qq and tetraquark interpolating 

^ . fields. Several difficulties which complicate the extraction of scalar meson masses from the 

■ ' scalar correlators are pointed out. One of the major difficulties is the presence of the two- 



pseudoscalar scattering states, which often dominate the correlator at light quark masses. 
In the simulations with unphysical approximations such scattering contributions are even 
more disturbing since they are often large and sometimes negative. Techniques which allow 
extraction of scalar meson masses in presence of scattering states are listed. Preliminary 
results of various lattice collaborations are presented. 

fi; 1 Introduction 



^ . The nature of light scalar mesons is not understood yet. Now there is a strong experimen- 

ON ! tal evidence for a light a meson, whose pole was extracted with a small uncertainty from 

experimental data using model independent Roy analysis [1]. Experimentally well estab- 
CN I lished scalar resonances below 1 GeV are also ao(980) and /o(980), while light fi;(800) is still 

■^ ■ controversial. At present it is still not clear whether these resonances represent the conven- 

tional qq nonet (its natural spectrum is illustrated in Fig. [T^ for vector case) or perhaps 
O ■ the tetraquark states. The diquark anti-diquark tetraquarks [q'g] 3^,3^ [^^] 3^,3^ are theoretically 

especially well motivated: they lead to inverted mass spectrum illustrated in Fig. [TJd, which 
has considerable similarity to observed spectrum [21 [3]. The physical states may in reality 



^ ■ be mixtures of qq and tetraquarks, while the singlets may mix also with glue. 



Lattice QCD has proven to be a very successful tool for determining the ground state 
hadron masses from first principles. The scalar meson puzzle would benefit from determina- 
tion of scalar meson spectrum on the lattice, but I will argue that this is more challenging 
than determination of masses for most of other ground states. The basic question that lattice 
simulations raise in this respect is what is the mass of the lightest state that couples to a 
given qq or tetraquark interpolator. The major problem is that the lightest scalar state is 
usually a two-pseudoscalar scattering state. The information on the one-particle scalar reso- 
nance comes from sub-leading contribution and its extraction is still challenging for present 
lattice simulations. 



^Talk presented at Scalar mesons and Related Topic (SCADRON70) in Lisbon, Portugal, February 2008. 



In this paper I will briefly review results from the simulations of light scalar mesons on 
the lattice. Two nice reviews on this topic were recently written by McNeile [1]. 

,- - uu/dd T - J KK 
miss - udss — ' < 
ss [_ Tir] 



Mass 



KK 



K* — > Kii usdd — . Kit 



10 > TTTTTT 



ud p — * TTTT udud 



1 -10 

1(a) 1(b) 



Figure 1: The illustration of the "natural" spectrum for (a) qq nonet, shown for the case of vector 
mesons, and (b) scalar tetraquark nonet. 



2 Methods for calculating scalar masses and related 
difficulties 

The basic object that is calculated on the lattice in order to extract the hadron mass is a 
hadron correlation function. Let me present the standard method for determining a hadron 
mass on the example of oq meson with 1 = 1. This example will also point out the specific 
difficulties to determine the scalar meson masses due to the presence of scattering states. 
Correlation function corresponds to creating a state (1(0, 0)m(0, 0) at t = and annihilating 
it at some later time t 

C{t) = ^((J(f,t)M(f,t)|M(0,0)rf(0,0)) , (1) 

X 

where total momentum is taken for concreteness. The interpolators with quantum numbers 
J^*" = O"^"*" and 1 = 1 may have various shapes as long as they couple with qq. For 
concreteness I have written the interpolators with point source and point sink. One of the 
main reasons why the scalar meson masses have not been yet reliably determined on the 
lattice is that the lattice correlators for light scalars and significantly more noisy than for 
light pseudoscalars and vectors. The second major difficulty has to do with the presence of 
open or nearly open decay channels: the created state d(0, 0)m(0, 0) with u and d at the same 
point does not correspond to pure oq meson. It creates all eigenstates of the Hamiltonian 
with quantum numbers J^*" = 0"'""'". In simulations of proper QCD with u, d, s dynamical 
quarks these are 

\d{0, 0)m(0, 0)) = A\ao) + B\a*o) + a\nr])i + Di\KK)i + Ei\7ir]')i + ... , (2) 



where the subscript i refers to various momenta of pseudoscalars. The time evolution of the 
correlator is obtained by inserting the full set of Hamiltonian eigenstates \n) 

C{t) = ^(|Jn|n)|2 e-^"* = |(Jn|ao)| V'"""* + ^ |(Jn|7rr/),| V^V + ... (3) 



with E^^ = m^r + rrirf, y ml. + kf + -w ?Ti^ + kf , .... Here ki = 27rj /L on the lattice with 
finite extent L. The scattering states with different momenta are well separated in energy 
for current lattice sizes. 

If the lightest energy eigenstate was Oq, then the extraction of its mass from C(t) would 
be simple: one would just have to fit niaQ from C(t) oc 6"™"°* at large t, where states 
with larger energies have exponentially decreased. In nature, however, ?72^ + ?7i^ is smaller 
than niao and the large-time behavior of the correlation function is dominated by the ttt] 
scattering state. The lightness of two-pseudoscalar scattering states therefore makes scalar 
meson spectroscopy particularly difficult. The scalar meson mass has to be extracted from 
the sub-leading contribution to the correlator. The available methods to do that include: 

• In the case of the scalar correlator with point source and point sink one can predict 
the two-pseudoscalar scattering state contribution using the corresponding version of 
Chiral Perturbation Theory (ChPT) up to a certain order of chiral expansion. The 
scalar meson mass can be determined by subtracting the predicted scattering contri- 
bution from the correlator and fitting that with e""*"''*. This method has been first 
applied in quenched QCD by Bardeen et al. [S]. Two fiavor and partially quenched 
QCD has been considered in [B] and three-fiavor staggered QCD in [7l[H]. The scalar 
simulations with chiral valence quarks on staggered sea were studied in [TJ [9l [10] , while 
twisted-mass QCD in [HI [12]. At light quark masses, the scattering states dominate 
the correlator and in all cases agree well with corresponding ChPT predictions. These 
are most often without any free-fit parameters, since they have all been pre-determined 
from other observables. 

• One- and two-particle states can be distinguished if correlator is computed at several 
lattice volumes since the spectral weight | (Jwln) p of a two-particle state is proportional 
to l/L^, while the spectral weight of one-particle state is not expected to be strongly 
volume dependent. Although the method is powerful, it has been only applied in [13] 
since results on various volumes are missing in dynamical simulations. 

• Scattering states can also be identified by their energy shift as a function of lattice size 
and pseudoscalar mass. The method is costly and was only applied in [21 [13] . 

• One can apply different boundary conditions for quark and antiquark fields, which 
elevate the threshold for two-particle state with respect to one-particle state 



A correlation matrix with various shapes of sources and sinks can be calculated instead 
of a single point-point correlator ([3]). In this way the energy of the ground and some 



of the excited states can be extracted using the variational method, used for scalar 
spectra by UKQCD [13 [ISl [IT] , BGR ^^ and RBC [2D] collaborations. 

In the following I discuss the scalar meson simulations which use various techniques described 
above. 

3 Simulations with qq interpolators 

ao meson 

First I consider the simulations with u, d, s dynamical quarks with relatively light 
u, d quarks, where one hopes to be simulating nature as closely as possible. The only ex- 
isting scalar meson simulations of this type use staggered valence quarks on a staggered sea 
(staggered QCD) [8] or chiral valence quarks on staggered sea quarks (mixed QCD) [9|. The 
mu4 masses in these simulations are so low that one expects the vrr/ scattering state to have 
the lowest energy, while the contribution of ao is expected to be sub-leading. Surprisingly, 
the ground state energies extracted from simulated correlators are not m,r + "^»?, but closer 
to 2?7i^ for both types of simulations when niu^d are sufficiently low [HI [9|. In addition, the 
correlator is negative in simulation that uses chiral quarks on staggered sea. 

Obviously, one is observing non-physical lattice artifacts, which get more and more pro- 
nounced as the u, d masses approach thier experimental value. In order to understand what 
is the source of these unphysical effects, the contribution of two-pseudoscalar intermediate 
state to / = 1 point-point correlators at finite lattice spacing was predicted in corresponding 
staggered and mixed ChPT [7]. It was found that these artifacts are due to the discretization 
effects, which arise via taste violation and the use of different actions for valence and see 
quarks. In the continuum limit all e"^'"'"* contributions were shown to cancel and analytic 
predictions agree with 7T'r],KK,7rri' scattering given by conventional ChPT [7]. At finite lat- 
tice spacing various e"^'"'"* contributions do not cancel and analytical predictions [7] agree 
well with simulated correlators in staggered [8| and mixed [9] QCD. Let me point out that 
these analytical predictions, which largely dominate correlators at low rriu^, have no free 
parameters at the lowest non-trivial order, since all of them have been determined by simu- 
lations of other quantities. The only free parameters of the fit are niao and spectral weight 
\{du\ao)\'^. The extracted masses are rriao = 1.01 ± 0.08 GeV [7] and niao = 0.87 ± 0.15 ^, 
but this can unfortunately not be very reliable values since the correlators are dominated by 
large and unphysical e"^'"'"* contributions. 

Even if the unphysical artifacts were not present, one expects that the correlators in 
simulations with u, d, s dynamical quarks would be dominated by e~(''"''+"^'')* at light rriu^d- 
For this reason I would propose that simulations with only u, d dynamical quarks have 
more chances to determine ruao, although they are less physical due to absence of dynamical 
s quarks. The only two-pseudoscalar intermediate state in this case is irr]', which is relatively 
heavy due to the anomaly and not so disturbing. The resulting masses from this type of 
simulations are shown in Table [H The key question is of course what is the mass of the 
lightest state that couples to du: does it correspond to ao(980) or ao(1450)? The first three 
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rriao 


SCALAR Coll. [21] 


u,d 


~ 2.4mp 


Prelovsek et al. [RBC] [6] 


u, d 


1.58 ± 0.34 GeV 


Prelovsek et al. [RBC] [B] 


u,d (part, quen.) 


1.51 ±0.19 GeV 


Hashimoto and Izubuchi RBC j20j 


u, d 


1.111 ±0.081 GeV 


McNeile, Michael, Hart 2003 UKQCD [I6] 


u, d 


1.0 ±0.2 GeV 


McNeile and Michael 2006 UKQCD [15] 


u, d 


1.01 ±0.04 GeV 


BGR 2007 [I9] 


u,d 


~ 1 GeV 


Bardeen et al. [5] 





1.34 ±0.09 GeV 


BGR 2008 [18] 





~ 1.4 GeV 


Mathur et al p^ 





1.42 ±0.13 GeV 



Table 1: Mass of the Ughtest 1=1 scalar meson from simulations with dynamical u, d quark and 
from quenched simulations. 



simulations [6l [2T| used a single correlator and were forced to fit at relatively low times, 
therefore the resulting mass might be artificially high due to the contribution of excited 
states. The rest of the two-flavor simulations [151 [HI [IHl 120] calculated a correlation matrix 
and determined niao using variational analysis, obtaining rriao — 1 GeV. The recent result 
from domain-wall fermions [20] is un update of [B] with larger statistics and 2x2 correlation 
matrix (using smeared and point source and sink): the resulting mass is lower, demonstrating 
the excited state contribution in the point-point correlator of [6]. McNeile and Michael [15] 
determined rrihi — rUao = 221(40) MeV on the lattice, compared to experimental result 245 
MeV, and deduced niao from experimental value of mbi. Although it appears that the more 
recent lattice simulations get a mass close to ao(980), the question whether ao(980) couples 
to qq is still open. The energies in the Table [1] simply correspond to the ground state energy, 
which may be Oq or vrr// - both have comparable energy. Before final conclusion on rriao 
may be drawn, the two flavor simulations have to identify ao as well as Tirj'. An attempt to 
estimate contribution of ttt]' to the point-point correlator was made in [6]: the ChPT renders 
the Trrj' contribution to be one order or more smaller than the lattice point-point correlators 
for ?72q ■^~ > 600 MeV in the fitted time range. 

One would naively think that the quenched correlator would be free of scattering states 
due to the lack of the sea quarks, but in this case the situation is even worse. The e"^"''"* 
contributions, which are large and negative, arise from unphysical irt]' scattering in quenched 
QCD [5]. In Table [1] I collect some of the recent quenched results for rriao which take these 
unphysical effects into account. A one-particle ao state was found in addition to two ghost 
scattering states in a variational analysis [22] of a 8 x 8 scalar correlator matrix in |18j . 
Analysis of a single point-point overlap correlator allowed the authors of [13] to extract two 
ghost states and a physical ao using a sequential empirical Bayes method. The quenched 
results for mao are higher than recent dynamical results. 

K and Kq mesons 

The relevant us correlator receives contribution from strange scalar meson as well as ttK 



scattering. The ttK state is so light that it has to be disentangled in any kind of simulation 
before conclusion can be drawn about tti^. So far the only simulation which demonstrated 
one and two particle states in us correlator was a quencheqj overlap simulation [13] , giving 
tuk* = 1.41 ±0.12 GeV. 

a meson 

Calculating the mass of the lightest singlet qq meson on the lattice is even more challeng- 
ing than for the nonsinglet mesons, described above. The vrvr scattering state is particularly 
light and has to be disentangled. The correlator is a sum of connected quark diagram, which 
is very sensitive to any non-physical approximations undertaken in simulations, and the dis- 
connected diagram, which is very expensive to evaluate numerically. On top of that the 
singlets mix with glueballs and their masses have shown sizable lattice spacing dependence 

A large scale simulation of a single qq correlator with u, d dynamical quarks on a fairly 
small lattice was performed by the SCALAR collaboration, giving 772^ — mp with a single 
exponential fit [21] • The extracted ttIq- is (slightly) below their vrvr threshold. In order to 
be confident that observed ground state corresponds to a, a further study is needed which 
would identify the vrvr scattering state. 

Lee and Weingarten [23] studied the mixing of qq and glue in a quenched simulation with 
quark masses around strange quark mass. They found that /o(1710) meson was 74(10)% 
glueball. More recently a dynamical simulation of qq and glue mixing was performed by 
UKQCD [15| using the variational analysis. They find the mass of approximately 1 GeV 
for the lightest isosinglet state, which couples to both qq and gluonic interpolators. In 
comparison, the quenched simulations render glueball masses well above 1 GeV. 

The singlet qq point-point correlator was simulated also with staggered valence quarks 
on a staggered sea [8]. As expected, the lattice correlator is largely dominated by vrvr scat- 
tering and agrees well with the tttt scattering prediction from staggered ChPT [S] . The only 
unconstrained parameters in the fit are related to sub-leading a contribution, leading to 
nia ~ 0.75 ± 0.16 GeV. The resulting rriu can unfortunately not very reliable due to the 
dominance of the vrvr scattering, whose contribution is even more complicated due to the 
taste breaking at finite lattice spacing. 

ETM collaboration simulated qq meson with two flavors of dynamical twisted-mass quarks 
and found a mass close to 2mjr [12], which is probably due to the tttt scattering. 

4 Simulations with tetraquark interpolators 

So far all lattice simulations with tetraquark interpolators have been quenched. They also 
all omitted the contribution of the disconnected quark diagrams to the correlator. These 
two approximation can be even justified in this case to some extent: one can unambiguously 
assign a quark content to a state by excluding processes that mix qq and (fq^ - dynamical 
quark loops and disconnected quark diagrams do exactly that. 



^Quenched / = 1/2 correlator receives e ('"t+™^)* due to the unphysical Krj' scattering state. 



Mathur and collaborators [13] calculated the 1 = correlator with interpolator (q'y^q) (q'y^q) 
at the source and sink. They used overlap quarks with small m^.d? corresponding to m^r down 
to 182 MeV. They analyzed correlators at spatial volumes 12^ and 16^. At each volume, they 
were able to extract three states from a single point-point correlator using the sequential em- 
pirical Bayes method. The ground state had energy slightly lower than 2771^^, while the ratio 
of the spectral weights for two volumes was consistent with 16^/12'^ (see Fig. [2]). This indi- 
cates a vrvr scattering state with zero relative momentum. The energy shift of this scattering 
state from the value 2r?7,^ agrees with the analytical prediction from quenched ChPT [21] at 
both volumes. The third state had an energy close to the energy of two pions with lowest 
non-zero momentum 27r/L. In the middle they found a state with energy well separated from 
any vrvr scattering state and with a ratio of the spectral weights for two volumes consistent 
with one. This gives an indication of a one-particle tetraquark state with a mass of about 
600 MeV, close to observed a [T] . This is one of the few available lattice studies which was 
able to separate the one and two particle contributions to the correlator. 
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Figure 2: Three lowest energy states extracted from a single quenched tetraquark correlator [13] 
are shown on the left. The ratio of the spectral weights for two volumes in case of lower two states 
are shown on the right. 



Alford and Jaffe [2] simulated the correlator with the same {q'j5q){q'y5q) interpolator in 
non-exotic 1 = and exotic 1 = 2 channels. The correlator was computed at a single 
relatively heavy quark mass and for a range of lattice sizes L. They studied the energy shift 
6E = E — 2?7i7r as a function of L, where E is the ground state energy extracted from the 
lattice correlator. They found that SE{L) for the exotic channel agrees with the fulo ChPT 
prediction for tttt scattering at finite volume. This indicates that the 1 = 2 ground state 
observed on the lattice is tttt scattering. The energy shift for 1 = does not follow the full 
ChPT analytic prediction for tttt scattering, providing another indication for the existence 
of a / = tetraquark. 

The Japanese group [H] made a simulation with diquark anti-diquark interpolator, which 

state in quark mass range rris < nig < 



is very well theoretically motivated [21 [3]. The / 



^It is surprising that the energy shift agrees with full ChPT prediction as one would have to compare it 
with quenched ChPT prediction. I guess this is due to the sizable mass of the quarks in [2]. 



2ms is simulated. They used conventional boundary conditions (periodic for q and q) as well 
as hybrid boundary conditions (anti-periodic for q and periodic for q). The second option 
does not alter tetraquark energy while it raises two-pion threshold (both pions have nonzero 
momenta) [H]. They extracted a ground state energy from a single correlator and it was 
consistent with tttt scattering in both cases. This indicates that a tetraquark resonance lies 
above the vrvr threshold, if it exists in the quark mass region m^ < mg < 2ms. 

5 Conclusions 

Lattice QCD has not yet given a definite answer what are the masses of the lightest qq and 
tetraquark resonances with J^ = 0"''+ and / = 0, 1/2, 1. The major reason why the consensus 
has not been reached is that two-pseudoscalar scattering states usually represent the state 
with the lowest energy in the simulations. The information about the resonant states has to 
be extracted from sub-leading contributions to the correlators, which is still challenging for 
most of the current lattice simulations. A quenched lattice simulation which was able to get 
an indication for a tetraquark a resonance in addition to tttt scattering state was presented 
in [13] (see Fig. [2]). The unquenched simulations of correlation functions with a variety of 
sources and sinks on a range of volumes would allow separation of scattering and resonant 
states, therefore allowing determination of spectra for ground state scalar resonances. I 
believe this is a feasible project for the near future. 
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